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ABSTRACT 

The technique of determining the acoustic location of layers of sharp changes in the sound speed inside a star from the oscillatory 
signal in its frequencies is applied on a solar-type star, the CoRoT target, HD 49933. We are able to determine the acoustic depth of 
the second helium ionisation zone of HD 49933 to be 794+gg s. The acoustic depth of the base of the convective zone is found to be 
1855!™ s where the large error bars reflect the ambiguity in the result, which is difficult to determine with present precision on the 
frequencies because of the intrinsically weak nature of the signal. The positions of both these layers are consistent with those in a 
representative stellar model of HD 49933. 

Key words. Stars: individual: HD 49933 - Stars: oscillations - Stars: interiors - Stars: fundamental parameters 



1. Introduction 

Any localised feature in the sound speed inside a star, such as 
that caused by the change in the temperature gradient at the base 
of the convection zone, introduces an oscillatory term in fre- 
quencies as a fu nction of radi al order n, that is proportional to 

1 1 pR fi ■ 

1990) where th = — is the acoustic 



sin(2rdW nj f + <f>), (Gough 



depth of the localised feature, c is the speed of sound, r& is the 
radial distance where the feature is located, io n { is the angular 
frequency of a mode with radial order n and degree (, and <p is 
a phase factor. This oscillatory signature has been extensively 
studied for the Sun in order to determine the extent of overshoot 
below th e solar convection zone ([Monteiro et al.| 1994] |Basu et| 
al. 1 1994 Roxburgh & Vorontsov 1994]. It has been proposed 
earlier that this may be used for distant stars also to find the 
position of the base of the convective envelope or the second he- 
lium i onisation zone (|Monteiro et al.|2 000 ; Mazumd ar & Antia] 
" 2001] |Roxburgh & Vorontsov||2003| |Ballot et al.|2004[|Basu et| 



al.| |2004| |Mazumdar p005| l. Indeed, |Miglio et al.| ( |2010| l have 



used the modulation of the frequency separations to determine 
the location of the second helium ionisation zone in a red giant 
star. 

We present here the results of applying this technique for 
the first time on a main sequence solar-type star: HD 49933, ob- 
served by the CoRoT mission in 2007-08 ( |Baglin et al.|2006| >. 

2. The technique 

The oscillatory signal in the frequencies is quite small and is 
embedded in the frequencies together with a smooth trend aris- 
ing from the regular variation of the sound speed in the stel- 
lar interior. It can be enhanced by using the second differences 
6 2 v(n, = v(n - \,€) - 2v(n, €) + v(n + 1, £), instead of the fre- 



quencies themselves (see, e.g., Gough 1990 Basu et aL] 1994 
Mazumdar & Antia|200l]|Basu et al!|2004| i. 

The acoustic depths of the base of the convective zone (BCZ) 
and the second helium ionisation zone (HIZ), tbcz and the, re- 
spectively, can be obtained from the data by fitting the second 
differences to a suitable function representing the oscillatory sig- 
nals from these layers (Mazumdar & Antia |2001| >. We follow 
Basu et al. ( 2004| l in choosing the following functional form: 



a Q + b sin(47rvT B cz+0BCz) 
+(c +c 2 /v 2 ) sin(47rvT HIZ +0Hiz), 



(1) 



where a Q , b , c , c 2 , t BC z, <Pbcz, ^hiz and H iz are 8 free param- 
eters of fitting. Since the number of data points in the observed 
set of frequencies for HD 49933 is quite small, we optimise the 
number of free parameters to strike a balance between a fair rep- 
resentation of the oscillatory signal and a reasonable^ 2 . We find 
that the smooth component is fairly constant over the range of 
frequencies that we use, as is the amplitude of the signal from 
the BCZ. However, the amplitude of the signal from the HIZ 
varies more sharply with frequency, and thus requires at least 
one frequency-dependent term. We note that Basu et al. ( 2004) l 
have shown that the exact form of the amplitudes of the oscil- 
latory signal does not affect the results significantly. To test our 
results, we also fit a function suggested by |Houdek & G ough 
( 2007|l (simplified version of Eq. (22) in that paper, with 8 free 
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parameters). 

The fit is carried out through a nonlinear x 1 minimisation, 
weighted by the errors in the data. The errors in the second dif- 
ferences are correlated, and this is taken care of by defining the 
X 1 using a covariance matrix. The effects of the errors are consid- 
ered by producing 1000 realisations of the data, where the mean 
values of the frequencies are skewed by random errors corre- 
sponding to a normal distribution. The successful convergence 
of such a non-linear fitting procedure is somewhat dependent on 
the choice of reasonable initial guesses. To remove the effect of 
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initial guesses affecting the final fitted parameters, we carry out 
the fit for multiple combinations of starting values. For each re- 
alisation, 100 random combinations of initial guesses are tried 
for fitting the function above and from among the resulting fits, 
the one which produces the minimum value of x 2 is accepted as 
the fit for that particular realisation. 

The median value of each parameter for 1000 realisations is 
taken as its fitted value. The error in the parameter is estimated 
from the range of values covering 68% area about the median 
(corresponding to lcr error). As an example, the histogram of 
distribution of the parameters tbcz and thiz over multiple real- 
isations are illustrated in Fig.[T] Thus the quoted errors in these 
parameters reflect the width of these histograms on two sides of 
the median value. 

The strength of the oscillatory signals are measured through 
the amplitudes of the respective oscillatory components in the 
frequencies themselves. These are denoted by A B cz an d ^hiz 
respectively, and estimated by averaging the amplitudes of the 
oscillatory components in Eq. [T] over the observed frequency 
range, and correcting for the scaling factor of 4 sin 2 (7rr/To) aris- 
ing in the second differences (see Mazumdar & Antia 2001 1. In 
the present study, we use the amp itudes only as an index of the 
significance of the detected depth of the acoustic glitches. 

We tested the technique on the frequencies of a polytropic 
model, with a polytropic index of 1.5, provided by |Christensen-| 
Dalsgaard & Mullan ( 1994) . Since there are no sharp acoustic 



Table 1. Results of fitting second differences of HD 49933 fre- 
quencies to Eq. [T] Median values of the acoustic depths and am- 
plitudes and their errors are listed. The last two rows give the 
calculated and the fitted values of a representative stellar model 
of HD 49933. 



features in such a model, we do not expect to find any oscillatory 
pattern in the second differences of the frequencies. Indeed, if we 
attempt to fit a function of the form given by Eq. [T]to the second 
differences of this model, we only find a signal corresponding to 
half of the large separation, which is the periodicity at which the 
successive points occur due to £ = or 2 and £ = 1. 

The range of acoustic depth that can be probed by this tech- 
nique is limited. The smallest depth that can be fitted corre- 
sponds to the range in the observed frequencies, and is ~ 350 s. 
On the other hand, the sampling frequency of the data points de- 
termines the largest acoustic depth that can be fitted. The average 
large separation of HD 49933 being 85 /t/Hz, the upper limit on 
the acoustic depth that can be fitted is ~ 2940 s. 

3. Results 

3.1. Fitting the CoRoT frequencies 

We use the second differences of the frequencies obtained by 
|Benomar et al.| ( |2009) > to extract the acoustic depths of BCZ and 
HIZ of HD 49933. We ignore frequencies which have errors of 
more than 3 jt/Hz. This gives a total of 3 1 second differences to 
be used for { = 0, 1,2. We consider four different subsets of the 
data, namely: (i) I = modes only (12 data points), (ii) ( = 1 
modes only (16 data points), (iii) t — 0, 1 modes (28 data points), 
and (iv) t — 0, 1, 2 modes (31 data points). This was done to test 
for consistency of the results when modes of different degrees 
are considered, especially since the t — and I — 2 modes 
in the data have, on the average, higher errors than the t — 1 
modes. However, the number of data points in the first two cases 
are only marginally larger than the number of free parameters 
in the fitting function, and these cases are used only for testing 
purposes. 

The fits to the second differences of the mean frequencies 
and the distribution of the fitted acoustic depths of BCZ and HIZ 
for different realisations of the data are shown in Fig[T] Although 
we fit 1000 realisations of the data, only those realisations with 
Tbcz ^ 2700 s are considered due to reasons explained below. 



t 


Tbcz (s) 


thiz (s) 


A BCZ (//Hz) 


A mz (//Hz) 





1877!™ 


g 47+ 248 
OH 1 -253 


1 -06-034 


, q 9 +5.80 


1 


1496!** 


812+^g 


75 +0 24 


-0.29 


0,1 


1626+262 


806+« 


0.66_ i 5 


1 40 +02S 


0,1,2 


1855+4], 


794t| 


71+019 

i -0.15 


i n.-)+0.34 
i - :,z -0.31 


Model 


1998 


882 






0,1,2 


1982+^ 


848!^ 


47 +019 

-0.13 


i Q7+0.24 
i ""-0.22 



A fit to the data of subset (iv) above to the form suggested by 
Houdek & Gough ( 2007 1 is shown also in Fig [T] Table [T] lists 
the median values of the acoustic depths with estimated errors 
and the corresponding amplitudes of BCZ and HIZ signals with 
respective estimated errors for the four data subsets. 

The acoustic depth of the HIZ obtained from the different 
sets are consistent within the errorbars. The amplitude of the HIZ 
oscillatory function is also significant at more than 2.5cr level in 
each set. 

However, the BCZ oscillatory signal is not so consistent 
among different sets. For all values of I, we get a peak in the his- 
togram close to the value to/2 ~ 2900 s, which corresponds to 
the Nyquist frequency. The occurrence of this peak could mean 
that the amplitude of the oscillatory function corresponding to 
BCZ is close to zero, so that the fitted sine wave has periodic- 
ity equal to the natural spacing of the data points, i.e., the large 
separation. While this seems to be the case for I = 1 modes, 
the second differences of the radial modes appear to have an 
actual oscillatory signal corresponding to this frequency! Thus 
this effect cannot be fully ascribed to the Nyquist effect. This 
is also borne out by the fact that this peak persists even when 
we combine data with multiple values of t, which should not be 
the case (see, e.g., Mazumdar & Antia 2001 ). We cannot find a 
satisfactory explanation as to the origin of this periodicity in the 
radial modes, apart from an artifact introduced due to the noise. 
In order to eliminate this artifact for all subsets, we discard the 
realisations for which the BCZ signal is fitted with an acoustic 
depth of around 2900 s. Few realisations yield tbcz values higher 
than the acoustic midpoint of the star (~ 2900 s), and these are 
discarded also. Actually, we consider only those realisations for 
which tbcz < 2700 s, since the peak in the histogram around the 
acoustic midpoint has a typical spread of 200 s on either side of 
the peak. This results in reduction of the number of valid reali- 
sations by nearly 30%. 

In all the cases, there is a peak in the histogram for tbcz 
around 1950s, which is most likely due to the base of the con- 
vective envelope in the star. For I = and I = 0,1,2, this is the 
dominant peak of the histogram, while for I — 1 and I = 0,1, 
this is the secondary peak. For the latter two cases, the dominant 
peak is centred around 1370 s. If this value is to correspond to 
tbcz it would imply a very shallow convection zone in the star, 
which is an unlikely scenario given the expected mass and age 
of HD 49933. This peak cannot be associated with any artifact 
due to the sampling frequency, but on the other hand, we cannot 
assign any physical significance to it. 

To understand the difference in the values of tbcz obtained 
for different t, we carried out the following test. For each case, 
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Fig. 1. Fits of Eq.[T|to the second differences of the mean CoRoT frequencies for HD 49933 and the histograms for the fitted values 
of the and tbcz for different realisations of the data are shown for different subsets (corresponding to different choices of ( values). 
The second differences of the frequencies of I — (blue), t = 1 (red) and 1-2 (green) modes of HD 49933 and their fit to Eq.fT 
(magenta curve) are shown in the left panel. The dotted cyan line in the bottom left panel shows the fit by the Houdek & Gough 



( |2007| l form. The corresponding histograms of the fitted values of tbcz (in red) and the (in blue) for different realisations are shown 
in the right panel. N denotes the total number of valid realisations after rejecting fits with tbcz > 2700 s. The solid bands at the top 
of the right panel indicate the range of initial guesses for the two parameters in each fit. 



we held the parameter tbcz in Eq. [T] fixed at a certain value, 
while allowing all other parameters to be free, and monitored 
the value of the reduced \ 2 obtained by the resulting fit to the 
mean second differences. The fixed value of tbcz was changed 
continuously between 1000 s and 4800 s to find the minima in 
the reduced x 1 as this parameter varies. The x 2 is plotted as a 
function of tbcz in Fig. [2] We find that the t = 1 modes show a 
deep minimum around 1370 s, but only a slight inflexion around 
1950s. In contrast, the location of the minimum and the inflex- 
ion are interchanged for the £ — modes. When these modes 
are combined along with the I — 2 modes, we have minima at 
both the values. The third minimum around 2900 s for all cases 
has already been discussed. The almost mirrored shapes of the 
graphs above tbcz values of 2900 s merely illustrate the aliasing 
effect, as elaborated in Mazumdar & Antia (200 1). 



3.2. Fitting model frequencies 
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Fig. 2. The variation of y 1 with the adopted value of tbcz is 
shown for four subsets of the CoRoT data on HD 49933 corre- 
sponding to different choices of the degree, I. The vertical dotted 
line shows the acoustic mid-point of the star. 



We fit Eq. [T] to the second differences of frequencies calcu- 
lated from a representative stellar model for the star HD 49933 
as well. The model is generated with the CESAM2k evolution 



code (Morel & Lebreton 2008 1 and uses standard physical in- 
puts, including convective core overshoot, but no diffusive mix- 
ing. While we use the theoretical frequencies as the mean values 
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Fig. 3. Fits of Eq. [T] to the second differences of the frequen- 
cies of a theoretical model for HD 49933 (top panel) and the 
histograms for different realisations of the data for the fitted val- 
ues of the and tbcz (bottom panel) are shown. See Fig. [T] for 
explanations of colours and symbols. 



in this case, we set the errors in the frequencies equal to the er- 
rors in the observed frequencies for the corresponding mode, in 
order to mimic the uncertainties present in the data. The results 
of such a fit are displayed in Fig. [3] The model values as well as 
the fitted median values of thiz and tbcz are given in Table [T] 

We find that the median values of the acoustic depths of both 
BCZ and HIZ are recovered from the frequencies. However, the 
error in the value of tbcz is quite large due to the fiat nature of 
the histogram in Fig. [3] In contrast, the histogram of the HIZ 
signal is sharply peaked and hence the estimated error in the me- 
dian value of thiz is also smaller. We conclude that the intrinsic 
oscillatory signal arising from the acoustic glitch at BCZ is quite 
weak, and it is difficult to extract this signal even with model fre- 
quencies. It is beyond the scope of this study to investigate why 
this signal is so weak, and it needs to be probed whether this 
is a common feature of all stars in the relevant part of the HR 
diagram. The HIZ signal, however, is strong and can easily be 
used to extract the acoustic depth of the second helium ionisa- 
tion zone. 



4. Summary 

We have applied the method of determining the acoustic depth of 
layers of sharp variation in sound speed in stellar interiors from 
the oscillatory signal in frequencies to the CoRoT primary target 
HD 49933. 

The oscillatory signal is made more pronounced in compar- 
ison to the smooth change in sound speed by using the second 
differences of the frequencies. We fit a function with two os- 
cillatory components corresponding to the two major layers of 
sharp variation in sound speed, viz., the base of the convective 
envelope and the second helium ionisation zone. To prove that 
the technique does not produce artificial signatures of oscilla- 



tory signals when there are none, we tested it on the frequencies 
of a polytropic model where there are no real layers of sharp 
variation in the sound speed. In this case the only significant 
oscillatory signal corresponds to the uniform spacing of the fre- 
quencies, as expected. 

We find that the oscillatory signal arising due to the sec- 
ond helium ionisation zone is quite strong in the CoRoT data 
of HD 49933. We estimate the acoustic depth of this layer to be 
794*5g s. This is borne out by the frequencies of modes of all 
available degrees: £ — 0, 1, 2, when considered either separately, 
or in combination. The acoustic depth of this layer in a represen- 
tative stellar model is consistent with this value, and can also be 
easily extracted by our technique. 

The oscillatory signal due to the base of the convective en- 
velope is far weaker in the data, in comparison to the helium 
signal. Moreover, we find that the signal is not consistent be- 
tween modes of different degrees. We find additional periodici- 
ties in the frequencies corresponding to different acoustic depths 
in modes of I — and € = 1, for which we are unable to as- 
sign any physical acoustic feature inside the star. It seems likely 
that these are artifacts due to the errors in the data which have 
no physical relevance. Nevertheless, we do find an oscillatory 
signal corresponding to an acoustic depth of 1855^2 s in the 
combined frequency set of different degrees, which is consistent 
with the position of the base of the convective envelope in a rep- 
resentative stellar model. This oscillatory signal is weak even 
in the model frequencies, and is fairly difficult to extract from 
the second differences. These results prove that the technique of 
extracting acoustic depths of sharp features inside a star can be 
successfully applied for real asteroseismic data. 
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